Adaptation of the Modified Adiabatic Approximation to Strong Field Ionization 



On 
O 
O 
(N 

£ 

1—5 



o . 

O 1 

CZJ ■ 

>v 

^ : 



> 

oo 

(N 

O 

On 
O 



X 



Denys I. Bondar, 1,2 'El Wing-Ki Liu, 1 '!]] and Gennady L. Yudin 2,3, [l 

1 University of Waterloo, Waterloo, Ontario N2L 3G1, Canada 
''National Research Council of Canada, Ottawa, Ontario K1A 0R6, Canada 
3 Universite de Sherbrooke, Sherbrooke, Quebec J1K 2R1, Canada 
(Dated: June 6, 2009) 

Savichev's modified adiabatic approximation [Sov. Phys. JETP 73, 803 (1991)] is used to obtain 
a general form of a quantum-mechanical amplitude of ionization by a low-frequency laser field. The 
method possesses only one requirement that the frequency of the laser field must be low. Connections 
of the obtained result with the quasi-classical approximation and other previous investigations are 
discussed. 
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Quantum-mechanical phenomena induced in a strong 
laser field are explicitly non-perturbative in nature 
Needless to say, theoretical investigations of such 
strong field phenomena go beyond the conventional 
perturbation theory. Beside numerical integration of 
the Schrodinger equation, many varieties of the cruasi- 
classical (QC) approximation are widely employed [l], [U, 
3]. In the current paper, we focus our attention on an- 
other approach - the adiabatic approximation. 

Roughly speaking, the adiabatic approximation can be 
introduced as follows. Once the frequency of the ex- 
ternal laser field is much lower than the characteristic 
atomic frequency, u> -C u> a t, an approximate solution of 
the Schrodinger equation can be found by means of aver- 
aging over atomic (internal) degrees of freedom. There- 
fore, the adiabatic approximation is a method of con- 
structing the asymptotic expansion of the solution of the 
non-stationary Schrodinger equation in terms of the small 
parameter u>/u> a t- 

Born and Fock [3| founded the theory of the adiabatic 
approximation for a discrete spectrum by formulating the 
adiabatic theorem. Landau [5| estimated the probability 
of nonadiabatic transitions between discreet states. How- 
ever, the leading-order asymptotic result for such a quan- 
tity was obtained by Dykhne Q. Afterwards, nonadia- 
batic transitions between discreet states was thoroughly 
analyzed by many authors (for reviews see Ref. [§]). 
Summarizing research in this area, one may conclude that 
nonadiabatic transitions in discreet spectra are quite well 
studied. 

Nevertheless, transitions from a discreet state to a con- 
tinuum one are the subject of on-going investigations for 
many decades. Despite much work in this topic, a univer- 
sally accepted approach is lacking. The direct general- 
ization of the Dykhne method was developed by Chaplik 
[Toj | . Exactly solvable models were reported by Demkov 
and Oserov [ll[, Ostrovskii [13], and Nikitin [l| (for 
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review see, e.g., Ref. [l4j). The advanced adiabatic ap- 
proach was introduced and widely employed by Solov'ev 
[la ]. Finally, Tolstikhin recently developed a promising 
version of the adiabatic approximation for the transitions 
to the continuum [l6| by employing the Siegert-state ex- 
pansion for nonstationary quantum systems [l7j . Yet, 
the approach presented in Ref. [161 ] is limited to finite 
range potentials. 

Our investigations are based on a seminal result that 
ought to be summarized foremost. Following the Solov'ev 
advanced adiabatic approach [ilj . Savichev [H[ proved 
the following. If the adiabatic state \tpi(ip)) and the cor- 
responding adiabatic term Ei(ip), 



H(<p)\ip i (<p))=E i {<p)\il> i (<p)) 



(1) 



are known, then the solution | of the nonstationary 
Schrodinger equation (the atomic units, h = m = \e\ = 1, 
are used throughout unless stated otherwise) 



(2) 



where ip = u>t is a phase of the laser field, subjected to 
the initial condition 



|^fe(B,)))e^^-<^[l + OH], 



has the following form within the adiabatic approxima- 
tion (u> -C 1) 



!*(*)> 



exp 




dEdtp' \i/>i(<Pi(E))) [l + 0(w)\ x 
<p' 



Ei{tp)(kp J - %Et 



(3) 



Here ifii(E) is the inverse function of Ei(ip). Note that no 
assumptions on a form of the Hamiltonian H were made. 

Let \i) and ]/) be stationary states (for specification 
see Eq. ([8]) and the comment after), and we shall assume 
that the quantum system with the Hamiltonian H is in 
the state \i) at t — — oo. The main aim of this section 
is to obtain the general form of the transition amplitude 
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9Jti_v/ that the given quantum system will be found in 
the state |/} at t = +00. 

Before going further, we are to introduce notations. 
First, we arbitrarily partition the Hamiltonian H: 



H(<p) = H (<p) + V(<p). 



(4) 



Second, we denote by /(*)) the solutions of Eq. © 
such that |*i(-oo)) = \i), |*/(+oo)) = |/); simi- 
larly, |3>»,/(£)) are the solutions of the nonstationary 
Schrodinger equation 

idt\$(t))=fi (<p)\*(t)), 

with the initial conditions: 00)) = \i) and 

|$/(+oo)) = I/}, correspondingly. 

Having defined all necessary functions, we introduce 
two equivalent forms of the transition amplitude 9Jtj_>/ 
by employing the corresponding version of the S-matrix 
(see, e.g., Ref. [19]): the reversed time form (sometimes 
called the "prior" form) 



= -i 



and the direct time form (the "post" form) 



(5) 



an 



(d) 



<$/(*)! y(wi)|*i(t)>dt. (6) 



It is noteworthy to recall the physical interpretation of 
Eqs. © and ©. The terms {$ f (t)\V(u>t)\$i(t)) and 
($/(i)| V{ojt) ys>i{t)) can be regarded as the amplitudes 
of quantum "jumps," which occur at the time moment t. 
The integrals over t convey that these jumps take place 
at any time. 

Introducing the adiabatic state \4>f(ip)) and term 
Ef(ip) of the Hamiltonian i? , 



& (<p)\M<p)) = E f(<p)\M<p)) 



(7) 



the wave function |$y(f)) can be readily presented in the 
form of Eq. ([3]). In further investigations, we employ the 
post form [Eq. ©], and thus we shall assume that 



|^(-oo)) = |t> , |^(+oo)) = |/) 



(8) 



In the case of the prior form [Eq. |(5J|], condition © has 
to be substituted by \4>i{— 00)) = \i) and |^/(+oo)) = 
I/), where \4>i{ip)) and \ipf((p)) are adiabatic states of the 
Hamiltonians Hq and H, correspondingly. 

Substituting the asymptotic representations [Eq. I|3])] 
of the wave functions [$/(*)) and [^(t)) into Eq. J6]), 
we obtain 



where z = (£,7],£',r)',<p) is a five-dimensional 
vector, <i 5 z = d£dr)d£'dr)'d<p, f(z) = 



(<l>f{tPf(s'))\v(tp)\M<Pi(£))h and s(z) 

s(v 7 <p) + s'b-rf) - rm)dn + r'Emdn. 

Bearing in mind that l/to is a large parameter, the 
five-dimensional integral in Eq. (J9j) can be calculated 
by means of the saddle-point approximation. Finally, 
the post form of the transition amplitude within the 
adiabatic approximation reads 



(d) 



2^ {M<p*)\v{<p*)\i>ii.<p*)) 



x exp 



[Ef(f) ~ Ei(ip)\ dkp\[l + 0(uj)}, (10) 



where denotes the summation over simple saddle 

points <p+, i.e., solutions of the equation 



±E f fa) * ±E M 



(11) 
(12) 



The physical interpretation of the sum over ip* is as 
follows: quantum jumps occur only at isolated time mo- 
ments t+ — tp+/uj, when the jumps are most probable; 
hence, are called "transition times." Note that the 
given interpretation deviates from the physical meaning 
of the time integral in Eq. ©. 

Some general remarks on Eq. (|10|) . the main result of 
this paper, are to be made: 

(i) (p* is usually a complex solution of Eq. (fTTj) ; there- 
fore, saddle points 92* with negative imaginary parts 
should be ignored because such points make exponen- 
tially large contributions to the amplitude, which leads 
to unphysical probabilities. 

(ii) If <pl, <p*,...,<p™ are solutions of Eq. (fTl]) , such 



that Im (ipi) > Im 



> . . . > Im (y>™) > 0, then all 



but the single term that corresponds to the saddle point 
92™ may be neglected in the sum over in Eq. (JTUJ) . One 
is eligible to do so since this saddle point has the largest 
contribution to the transition amplitude. 

(iii) On the one hand, the explicit form of Ef(ip) is 
solely determined by partitioning [Eq. Q]; on the other 
hand, Ei(ip) is unique for a given quantum system. 

(iv) The exponential factor of Eq. pH|) is similar to the 
exponential factor in the Dykhne approach [(| 0, El (see 
also Refs. [13, Hll) ~ the methods for calculating the am- 
plitude of bound-bound transitions within the adiabatic 
approximation. Hence, Eq. (|10p may be considered as 
a generalization of the Dykhne formula for bound-free 
transitions. 

(v) By employing an appropriate version of the saddle- 
point method, one can in principle generalize Eq. (JTDJ) 
for the case when condition ([12")) is violated. 

Now, the connection between the amplitude [Eq. (flU)) ] 
and the method of complex trajectories is to be man- 
ifested. According to the method of complex classical 
trajectories (see, e.g., Refs. HFj, [22|, and the imagi- 
nary time method [3(), to calculate the probability of the 
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transition from the initial state to the final, one should 
first solve the corresponding classical equations of mo- 
tion and find the "path" of such a transition. However, 
this path is complex; in particular, the transition point 
r* and transition time at which the transition occurs 
are complex. Parameters r* and are determined by 
the classical conservation laws. Next, one has to obtain 
the classical action Sf(rf,tf,r+,t+) + Si(r+, t+; r i: ti) for 
the motion of the system in the initial state from the ini- 
tial position Ti at time to the transition point at 
time and then in the final state from at t* to the 
final position Tf at time tf. Finally the probability of 
the transition is given by 

T tx exp{-2Im [S f (r f , tf, r*, t+) + <Si(r*, i*; r i( ti)]} . 

(13) 

Equations (|10|) and (| 13[) must coincide in some region 
of parameters. The method of complex trajectories can 
be derived as the QC approximation of the transition 
amplitude [Eq. ([5]) or Eq. l[6])] ; we outline this derivation 
below. Therefore, it would be of methodological interest 
to establish an explicit connection between Eqs. (fT0|) and 
(fT3|) . This connection has not been demonstrated. 

Without loss of generality, assuming that V(ivt) is a 
non-differential operator, we obtain the QC approxima- 
tion to Eq. © 



SOT 



(d) 



i dt d 3 vd 3 v f d 3 r l (f\ T f )V(ut, r) <r ( | i) 



xF f Fi exp {i [Sf(Tf,t f ; r, t) + S,(r, t; n, U)]} , (14) 

where tfj — ±oo, Ff exp(iSf), and Fiexp(iSi) are the 
QC versions of the propagators with the Hamiltonian Hq 
and H, correspondingly. We recall that the general form 
of the QC propagator is given by 



F {a) (r, t; r', t 1 ) exp \iS (a) (r, t; r', t') 



(15) 



where the sum denotes the summation over classical 
paths that connect the initial (r',t') and final (r, t) 
points. Therefore, usage of this form of the QC prop- 
agator, F exp(iS), is justified if we assume that there is 
only one such path; indeed, this is the case in the ma- 
jority of practical calculations, and thus we shall accept 
this assumption hereinafter. 

In order to reach Eq. f|13|) from Eq. (fT4"|) , one has to 
calculate the integrals over r and t in Eq. (fl"4|) by means 
of the saddle-point approximation. The equations for the 
saddle points and i.e., the transition points, read 

dtiS^t-^^U) ~ S f (r,t;r f ,t f )]\ t=ttr ^ = 0, (16) 

VrlSiir^r^ti) - S f (r,t;r f ,t f )]\ t=U r=it = 0. (17) 

Recalling the Hamilton- Jacobi equation 

d t Sij(v, t; Vij^ij) = -Wi,/(r, Vij,t), (18) 

where Hij are classical Hamiltonians and T'ij are clas- 
sical canonical momenta 

Vij(r,t) = V r Sij(r,t;rij,Uj), 



we rewrite Eqs. (|16l) and (|17|) as the law of conservation 
of canonical momentum and the law of conservation of 
energy: 



(19) 



W/(r*,7» / (r*,t*),^)=7< i (r*,7» i (r*,^),**).(20) 

Having introduced all the necessary quantities, we 
demonstrate the correspondence between Eqs. (fTU)) and 
(|13p within exponential accuracy. Performing a simple 
transformation and using Eq. (I18j) . we reach 

^/(r/.f/jr*,^) + Si(iv,^;ri,ij) 
= Sf(rf,tf,r*,ti) + Si(r*,ti;ri,ti) 

[d T Si(r±,T; r h ti) - d T S f (r±,T; r f ,t f )} dr 



[H f (r* , 7>/ (r* , r), r) - Hi (r* , (r* , r), r)] dr 



+Sf(rf,tf,r it) t i ) + 5i(r*,^;ri,ii 



(21) 



Usually in the case of multiphoton ionization, t* is com- 
plex and r* is real. Therefore, the last two terms affect 
only the phase and does not contribute to the probability. 

Finally, since TCf and Hi are the QC limits of Ef and 
Ei (this will be demonstrated below), we conclude that 
the exponential factors of Eqs. (fT0|) and (fl3|) indeed co- 
incide within the QC approximation. 

The wave function 

* gc (r,<)= / F i (r,t ] r',U)e^ Si ^ r '^cl>i n (r')d 3 T',(22) 



is the (leading-order term) QC solution of Eq. ([2]) with 
the initial condition $> gc (r,ti) = 4>i n (r). Employing Eq. 
(HHJ) and bearing in mind that Hi — Hi(r, T'i,t) does not 
depend on r', we obtain 



H^ qc = im^ 



[1 + 0(h)] = H^ qc [1 + 0(H)] . (23) 



Since we have freedom of choosing the initial condition 
<f>i n (r), there are in general infinitely many wave func- 
tions [Eq. ([221)] that satisfy Eq. (gSJ). Comparing Eqs. 
(|23| and (JTJ) , and taking into account the latter, we for- 
mulate the following property of the adiabatic term and 
state of a given quantum system within the QC limit: 
there exists only one adiabatic term, which is equal to 
the classical Hamiltonian, and any solution of the corre- 
sponding Schrodinger equation is also an adiabatic state 
that corresponds to this adiabatic term (i.e., the adiabatic 
term is infinitely degenerate). Note that this property is 
completely ruled out once the general form of the QC 
propagator (|15p is considered. 

The property stated above, nevertheless, merely accen- 
tuates the fundamental difference between the QC and 
adiabatic approximations. As mentioned in the intro- 
duction, the adiabatic approximation allows us to obtain 
the solution of the non-stationary Schrodinger equation 
as an asymptotic series in terms of the small parameter 
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Lu/uj a t; however, the QC approximation is a method of 
obtaining an asymptotic expansion of the solution with 
respect to the small parameter h. These two series are 
dissimilar in the general case. 

In conclusion, we clarify a connection between Eq. (|10[) 
and previous studies. Let H be the Hamiltonian of an 
atom in a low-frequency laser field, 

H(<p) = p 2 /2 - Z/r + V L (<p), (24) 

where the term Vz,(<p) describes the interaction between 
the laser field and the atom. In Ref. [2l[, the Keldysh 
theory [23] has been formulated within the Dykhne ap- 
proach by choosing Ef(<p) = A 2 (</j)/2 as the energy of 
classical oscillations of a free electron in the laser field 
and Ei(ip) = —I p , where A(ip) = — (F/w) sin<^ is the 
vector potential of the linearly polarized laser field, F 
is the strength of the laser field, and I p is the ioniza- 
tion potential. To obtain a general expression for single- 
electron spectra (within exponential accuracy), this for- 
mulation has been generalized in Ref. [24] by setting 
E f (ip) = [k + A(<^)] 2 /2 and E^ip) = -I p , where k 
being the momentum of the electron at the detector. 
Equation (jTDJ) rigorously justifies such employment of the 
Dykhne theory. Indeed, partitioning Hamiltonian (1241) 
such that V = —Z/r [see Eq. (J2J], we conclude that 
according to Eq. ((7]), \<f>j (ip)) is the Volkov wave func- 
tion; thus, Ef(ip) = [k + A(<ys)] 2 /2 is an exact equality 



within the given partitioning of the Hamiltonian. How- 
ever, Ei(tp) — —Ip can be obtained only when the dy- 
namical Stark shift is ignored. 

Let the Coulomb potential energy be represented as 

- Z/r = V shr {r) + V lng {r), (25) 

where the potential Vi ng (r) has a long-range behavior 
identical to —Z/r, but no singularity at the origin, and 
Vshrir) is a singular but short-range potential. Parti- 
tioning [Eq. ([25]) ] has been originally introduced in Ref. 
[25l | : recently, the similar partitioning has been exploited 
in the case of two-electron systems [26| . For the following 
representation of Hamiltonian |4]) : 

A (<p) = p 2 /2 + V lng (r) + V(ip) = V shr (r), 

Equation (10) is an adiabatic formulation of Eq. (23) of 
Ref. [25[. The last equation represents the amplitude of 
single-electron ionization that accounts for the Coulomb 
corrections. Therefore, Eq. (JTDJ) allows for further inves- 
tigations (beyond the QC approximation) of the impor- 
tant problem of the influence of a Coulomb potential on 
ionization rates @, 0, IM E3 ■ 
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